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The  rapid  growth  of  computer  communication  has  motivated  an 


intense  interest  in  packet  switching  radio  techniques  [1]. 
Furthermore,  there  is  a  growing  need  for  computer  communication 
and  information  distribution  in  tactical  military  applications, 
where  spread  spectrum  waveforms  must  be  used  in  order  to  achieve 
reliable  operation  in  the  precense  of  intentional  interference 
(jamming) .  As  a  result,  a  thorough  investigation  of  spread 
spectrum  packet  radio  networks  is  necessary. 

The  bit  error  probability  induced  in  frequency  hopped  spread 
spectrum  systems  has  been  examined  before  [2].  In  [2]  the  bit 
error  probability  is  computed  for  two  different  models  of  random 
frequency  hopping  patterns.  In  the  first  model  it  is  assumed  that 
each  random  frequency  hopping  pattern  is  a  sequence  of  independent 
random  variables  (i.e.  memoryless  frequency  hopping  patterns), 
while  in  the  second  model  it  is  assumed  that  each  random  frequency 
hopping  pattern  is  a  first  order  Markov  chain  (i.e.  Markov 
frequency  hopping  patterns) . 

The  computation  of  the  packet  error  probability  induced  in 
frequency  hopped  spread  spectrum  systems,  which  utilize  memoryless 
frequency  hopping  patterns  has  been  examined  before  ([3], [4], [5]). 
In  this  paper  we  are  going  to  compute  the  packet  error  probability 
induced  in  spread  spectrum  systems,  which  utilize  Markov  frequency 
hopping  patterns.  What  makes  the  problem  difficult  is  that  the  bit 
errors  are  not  independent.  Hence,  we  cannot  extend  the  results  in 
[2],  in  a  trivial  way,  in  order  to  compute  the  packet  error 


probability.  Nevertheless,  some  comparisons  with  the  packet  error 
probability  induced  if  we  assume  that  the  bit  errors  are 
independent  are  going  to  be  made. 

The  organization  of  the  paper  is  as  follows.  In  section  2  we 
describe  the  model  of  our  spread  spectrum  system.  Then,  in  section 
3  we  present  a  method  to  compute  the  packet  error  probability 
induced  in  our  system.  In  section  4  we  utilize  the  numerical 
results  of  section  3  and  an  educated  conjecture  to  evaluate  the 
performance  of  our  spread  specrtum  system.  Finally,  in  section  5 
we  give  a  brief  summary  of  the  accomplishments  of  the  paper,  and 
we  discuss  some  extensions  of  the  model  presented  in  section  2. 

2 .  The  Model 

The  frequency  hopping  system  will  now  be  described.  The 
channel  time  is  divided  into  slots,  and  the  users  in  the  network 
initiate  their  packet  transmissions  at  the  beginnings  of  slots. 

The  frequency  spectrum  is  divided  into  q  frequency  bins  and  the 
packets  are  divided  into  M  bytes  each.  Every  user  in  the  network 
sends  each  of  the  M  bytes  of  his  packet  at  a  frequency  bin  ,  which 
is  different  from  the  frequency  bin  used  by  the  previous  byte,  but 
equally  likely  to  be  any  one  of  the  remaining  q-1  frequency  bins 
(Markov  frequency  hopping  patterns) .  Furthermore,  different  users 
in  the  network  have  statistically  independent  frequency  hopping 
patterns.  We  also  assume  that  a  packet  consists  of  exactly  one 
codeword  from  a  Reed-Solomon  (RS)  code  for  which  up  to  e  byte 
errors  can  be  corrected.  A  packet  is  declared  successfully 


2 


transmitted  if  at  most  e  byte  errors  occur. 


3 .  A  method  to  compute  the  packet  error  probabilit's 


Let  us  assume  that  K  (K^2)  packets  are  transmitted  in  the 


same  slot.  These  packets  correspond  to  K  different  users  in  the 
network.  We  assign  indices  to  these  packets  (i.e.  packet  #  l, 
packet  #2,  ....  packet  #  K) .  Let  us  also  assume  that  the  receiver 


locks  on  to  packet  ft  1.  We  say  that  the  jth  byte  of  packet  #  1  is 
hit  if,  during  its  reception  by  the  receiver,  at  least  one  of  the 

other  packets  (i.e.  packets  #  2,  #  3,  . ,  #  K)  occupies  the 

same  frequency  bin  that  packet  #  1  occupies.  Let  us  now  denote 


by  <Pe(K) ,  the  probability  that  packet  #  1  is  decoded  incorrectly 
by  the  receiver  given  that  K-l  other  packets  interfere  with 
packet  #  1.  Our  objective  is  to  provide  a  method  for  the 
computation  of  pe(K) . 

We  denote  by,  {  f*  ;  l<j^M  },  the  frequency  hopping  pattern 
corresponding  to  packet  #  i  (l<i<K).  In  figure  1  we  show  a 
realization  of  the  K  packet  arrivals  at  the  receiver  site.  It  is 
worth  noting  that  the  realization  of  packet  arrivals  in  figure  1 
corresponds  to  the  worst  possible  case;  in  other  words  Pe(K)  is 
maximized  when  the  realization  of  figure  1  occurs.  This  is  true, 
because  with  the  realization  of  figure  1  ,  during  the  reception  of 
every  byte  of  packet  #  1  ,  all  other  K-l  interfering  packets  are 
present.  We  will  compute  Pe(K)  for  the  realization  of  packet 
arrivals  depicted  in  figure  1.  Furthermore,  we  will  make  the 
pessimistic  assumption,  as  in  [3],  that  when  a  byte  is  hit  a  byte 


error  results. 


We  will  compute  first  the  packet  error  probability  ,  pe(K), 
for  K=2 .  In  doing  so,  we  will  be  able  to  describe  the  major  points 
of  the  methodology  better.  Then,  we  will  discuss  the  steps 
required  to  compute  Pe(K)  for  Ks.3 . 

Case  1  .  K=2. 

It  is  obvious  from  (1)  that,  in  order  to  find  p  (2),  we  must 
compute  the  Pr[S(l,M)=i]  for  i=e+l , e+2 , . . . , M  .  From  the  formula  of 
total  probability  we  get: 


Proposition  1.  For  every  M£l  and  i^M  the  conditional  probability, 

112  2 

Pr (S (1 ,M) =i/f 1=s1, f 1=s1) r  depends  only  on  the  following  two  events 

.,12  .  . ,  1 ,  2 

i)s1=s1  or  n)s1+s1 


In  other  words  Pr(S (l,M)=i/f1=s1, f1=s1)  does  not  depend  on 

1  2  12  12 

the  actual  values  of  s1  and  s^  as  far  as,  s1=s1  or  s^s^. 

Proposition  1  is  proven  in  Appendix  A. 

Let  us  now  define 

fj=(fj,fj)  l<j<M  (3) 

The  proof  of  Proposition  1  is  based  on  the  following  Lemma 
Lemma  1 .  The  sequence  {  fj  ?  l<j<M  }  is  a  Markov  chain  with 
stationary  transition  probabilities. 

Lemma  1  is  also  proven  in  Appendix  A.  Due  to  proposition  1 
equation  (2)  becomes. 

Pr [ S ( 1 , M) =i ] =  q-1  Pr (S (1,M) =i/f|=l, f*-l) + 

+(l-q-1)  Pr (S ( 1 ,M) =i/f |=1 , f ^=2 )  (4) 

A  byproduct  of  the  proof  of  proposition  1  is  that  the 

.  .  .  .  .  1  2 
conditional  probabilities,  Pr  (S  (1,M)  =i/f  1=1,  f-L=l)  and 

1  2 

Pr  (S  (1  ,M)  =i/f  ^1 ,  f  ^2 )  ,  satisfy  certain  recurrent  expressions 
(see  Appendix  A  page  16  for  more  details) .  We  write  these 
recurrent  expressions  in  the  sequel. 


Pr(S(l,l)=0/fJ=l,f*=l)=0  (5a) 

Pr(S(l,l)=l/f*=l,f*=l)=l  (5b) 

Pr(S(l,l)=0/fJ=l,f*=2)=(q-2)/(q-l)  (5c) 

Pr(S(l,l)=l/fJ=l,f*=2)=l/(q-l)  (5d) 

Pr(S(l,n)=0/fJ=l,f*=l)=0  (6a) 

;2^n^M 


Pr(S(l,n)=i/f1=l,f‘=l)=(q-l)  *  Pr (S ( 1, n-1) =i-l/f ^=1 , f ‘=1) + 

;2^n<M,  0<i<n 

+(q-2) (q-1)”1  Pr(S(l,n-l)=i-l/f*=l,f*=2)  (6b) 

Pr(S(l,n)=0/f]=l,f*=2)=(q-2)2(q-l)-2  Pr  (S  ( 1 ,  n-1)  =0/f  ]=1 ,  f  \=2 )  (7a) 

;2<n<M 

Pr (S ( 1 , n) =i/f |=1 , f  ^=2 )  =  (q-2 )  (q-1)"2  Pr (S ( 1 , n-l)  =  i/f J=1 , f *=1)  + 
;2^n£M,  0<i<n 

+(q-l)-1  Pr (S (1, n-1) =i-l/ f t f ) + 

+ ( q~2 ) 2 ( q-1 ) ~2  Pr (S ( 1 , n-1) =i/f ^=1 , f 2=2 )  (7b) 

Pr (S ( 1 , n) =n/ f ^=1 , f ^=2 ) = (q-1) -1  Pr (S ( 1 , n-1) =n-l/f \=1 , f 2=2 )  (7c) 

;2^n^M 

Based  on  expressions  (5)  through  (7)  we  can  compute 
Pr (S ( 1 , M) =i/f |=1 , f ^=1)  and  Pr (S (1,M) =i/f |=1 , f2=2)  for  arbitrary  M 
and  i^M  as  follows. 

For  M=1  all  we  need  are  expressions  (5) .  For  M>1  we  start 

1  2 

from  n=2  and  we  evaluate  the  probabilities  Pr (S ( 1 , n) =i/f 1=1 , f 1=2 ) 

1  2 

and  Pr (S (1 , n) =i/f 1=1 , f ^=2 )  for  i=0,l,...,n  based  on  expressions 
(6)  and  (7).  Then,  we  perform  similar  computations  for  n=3 , 4 , . . . 
up  to  n=M.  Finally  we  end  up  having  computed  the  probabilities 
Pr(S(l,M)=i/f]=l,f2=l)  and  Pr (S (1,M) =i/f J=l, f2=2)  for 
i=0,l,2, . M. 

i  2 

Once  we  have  computed  Pr(S(l,M)=i/f  =l,f  =1)  and 
Pr (S (1 ,M) =i/f ^=1 , f2=2 )  for  i=0 , 1 , 2 , . . . , M  ,we  can  find  Pr(S(l,M)=i) 
for  i=0 ,1,2, ... ,M  through  formula  (4).  As  a  result  we  can  compute 
p  (2)  via  formula  (1). 


The  same  methodology  applies  when  K^3.  Hence, 


q  q 


Pr[S(l,M)=i]= 


Pr(S(l,M)=i/fJ=sJ1 


^K_  K 
,fl  Sl) 


-K 


(8) 


Furthermore,  a  proposition  and  a  Lemma,  which  are  similar  to 
proposition  1  and  Lemma  1,  can  also  be  stated  for  every  K>3 .  What 

makes  the  method  more  complicated,  as  K  increases,  is  that  the 

11  K  K 

number  of  events,  on  which  Pr  (S  (1,M)  =i/f  ^s^,  .  . .  ,  f  1=s1)  depends 
on,  increases  as  well.  For  example,  when  K=3  the  probability 
Pr (S (1 ,M) =i/f^=s|, f^=s^, f^=s^)  depends  on  the  following  four 
events: 


.  .  .  1  2  .  1  ,  3  1.2  .  1  3 

or  n)s1=s1  and  s1^s1  or  s^s.^  and  s1=s1 

12  13  2  3 

or  iii)s1^s1  and  s1^s1  and  s^s^^ 

12  13  2/3 

or  ivjs^s^^  and  s1^s1  and  s1rs1 

Nevertheless,  once  the  distinct  events,  on  which  the 

11  K  K 

probability  Pr (S (1,M) =i/f1=s1, . . . , f1=s1)  depends  on,  have  been 
correctly  identified,  recursive  expressions,  similar  to 
expressions  (5)  through  (7),  can  also  be  written  for  every  Ki3. 
Consequently,  following  the  same  methodology  presented  in  case  1, 
one  can  compute  Pe(K)  for  K^3 . 

Numerical  Results. 

In  Table  1  we  have  included  pg(2)  and  Pe(3)  for  q=10,25  and 


50  and  for  the  (31, 7) , (31 , 15) , (31 , 23) , (63 , 15) , (63 , 31)  and  (63,47) 


7 


Reed  Solomon  (RS)  codes.  In  the  same  table,  we  have  also  included 


§>  (2)  and  §>  (3) ,  where  p  (K),Ki2,  corresponds  to  the  packet  error 
probability  induced  by  our  spread  spectrum  system  under  the 
assumption  that  the  byte  errors  of  packet  #  1  are  independent. 

From  table  1  we  can  make  the  following  observations. 

1)  ?e(2)  is  an  upper  bound  of  Pe(2)  and  Pe(3)  is  an  upper 
bound  of  Pe(3)  for  most  entries  of  table  1. 

2)  iPe(2)-£e(2)|  /£e(2)  and  IPe(3)-gre(3)| /gie(3)  are  smaller 
than  one  for  all  entries  of  table  1. 

We  will  see,  in  the  next  section,  how  observation  2  will  be 
helpful  in  evaluating  the  performance  of  our  spread  specrum 
system. 

Note  that  if  K  users  are  present  in  a  slot  then, 

M 

Pe(K)=  ^  (”)  CPb(K)]1  tl-pb(K)]M_1  (10) 

i=e+l 

where  Pb(K)  is  the  byte  error  probability  given  K.  It  is  easy  to 
show  that 

pfa(K)=  l-(l-2/q)K"1  (11) 

where  q  is  the  number  of  frequency  bins  in  the  frequency  spectrum. 

4.  Performance  evaluation  of  the  spread  spectrum  system. 

Suppose  that  the  input  packet  arival  process  per  slot  is 
Poisson  with  intensity  s.  Then,  let  us  define  (as  in  [4]) 


to  be  the  average  packet  error  probability  induced  by  our  spread 
spectrum  network  at  input  rate  s  (note  that  Pe(0)=pe(l)=0) .  We 
also  define  (as  in  [4])  the  maximum  average  interference  level 
that  can  be  accomodated  at  a  given  packet  error  probability  v  by 

s*(v)=  max{ s :  qe(s)$v}  (13) 

★ 

Clearly  the  value  of  s  (v)  depends  on  the  number  of  frequency 
slots  q  and  the  code  rate  r=k/M  (k=number  of  information  bytes, 
M=total  number  of  bytes  in  the  code) .  If  we  account  for  the 
expansion  of  bandwidth  induced  by  q  and  r,  then  we  can  define  the 
normalized  maximum  average  interference  level  as 

S*(v)=  r  s*(v)/q  (14) 

Let  us  also  define, 

CP 

e~S  sK/K!  ?e(K)  (15) 

K=0 

(note  that  Pe(0)=Pe(l)=0) 

§f*(v)=  max{ s :  3e(s)<;v}  (16) 

and  £>*(v)  =  r  ^(vj/q  (17) 

According  to  observation  2  of  the  previous  section  the  ratios 
IPe(2)-P'e(2)  I  /Pe(2)  and  IPe  ( 3 )  -pe  ( 3 )  |  /Pe  (3 )  are  smaller  than  one. 
We  can  now  state  a  conjecture. 

Conjecture  :  |  Pe  (K)  -p  (K)  |  /f!  (K)  £  1  for  K£4  . 

If  this  conjecture  is  true,  and  we  find  s  such  that 


then  the  same  s  guarantees  that 

qe(s)  <  2v  (19) 

The  above  discussion  indicates  that  we  can  use  ?e(s)  and 
3*(v)  as  reasonable  estimates  of  the  performance  of  our  spread 
spectrum  system,  provided  that  the  above  conjecture  is  true.  It  is 
worth  noting  that  cJe(s)  and  3  (v)  are  easily  computable. 

5.  Conclusions. 

We  have  described  (section  3)  a  method  of  computing  the 
packet  error  probability,  Pe(K) ,  induced  in  a  spread  spectrum 
system  which  utilizes  first  order  Markov  frequency  hopping 
patterns.  We  have  also  computed  pg(K)  for  K=2  and  K=3 .  Then,  we 
made  some  comparisons  between  Pe(2)  and  ?>e(2)  and  between  Pe(3) 
and  pe(3) ?  pe(K),  K^2,  was  defined  to  be  the  packet  error 
probability  induced  by  our  spread  spectrum  system  if  we  assume 
that  the  byte  errors  are  independent.  Based  on  these  comparisons 
and  an  educated  conjecture  we  have  evaluated  the  performance  of 
our  spread  spectrum  system  (section  4). 

The  consideration  of  a  slotted  channel  in  section  2  is  not  so 
restrictive.  The  methodology  of  section  3  and  the  discussion  of 
section  4  are  still  valid  for  the  unslotted  channel,  provided  that 
we  resort  to  an  upper  bound  for  the  packet  error  probability. 
Consider  the  transmission  of  a  given  packet  in  the  unslotted 
channel,  and  consider  a  second  situation  in  which  it  is  assumed 
that  the  interference  level  is  constant  and  equal  to  the  maximum 


number  K*  of  interfering  transmissions  that  take  place  at  any  time 
during  the  transmission  of  the  given  packet  in  the  unslotted 
system.  Clearly,  as  table  1  also  indicates,  the  packet  error 
probability  for  the  second  system  is  larger  than  for  the  first. 
Hence,  pg(K*)  of  table  1  will  be  an  upper  bound  of  the  packet 
error  probability  of  the  unslotted  system. 

Furthermore,  the  pessimistic  assumption  that  a  byte  hit 
results  in  a  byte  error  need  not  be  made  either.  More  optimistic 
assumptions  described  in  [5]  section  IV,  where  thermal  noise  is 
also  present,  can  be  incorporated  in  our  model  too.  They  will 
simply  make  the  presentation  of  section  3  more  complicated. 
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Appendix  A 

We  will  first  prove  Lemma  1.  Let  us  now  denote 

;  j>l  (A. 1) 
as  in  (3) .  Furthermore,  we  denote 

•  j>l 

We  will  show  first  that 


1  2 

^j=(sj'sj) 


(A. 2) 


Pr(fj=sj/f._1=sj_1,...,f1=s1)=Pr(f.=sj/fj_1=sj_1)  (A. 3) 

for  any  s^,...,Sj  belonging  in  the  state  space  I  (in  our  case 

1=  {  (i1,i2):  i1=l,...q,  i2=l,...q  })  and  for  any  j^2. 

(A. 3)  is  equivalent  to 


Pr ( f j  =s j , . . . , f1=s1)=Pr(f j_1=Sj_1, . . . , f 1=s1 )  Pr(f j=£j , f 

cpr- (f  r1  (a. 4) 


But, 


Pr(fj=Sj, . . . , f 1=s1) =Pr (f l=s ) , f j=sj , . . . , f J=s|, f \=*\)p 

i  tt  2 

independence  of  { f ^ ; j > 1 }  and  {fjjj^l} 

11  11  22  22 

=Pr ( f .=s ! , . . . , f ’=s!)  Pr(ff=sT, . . .  ,f*=s^)  = 

=Pr (f 1=Sj/f j_1=s1_1, . . .  ,fj=sj)  Pr(fj_1=Sj_1, . . . ,f|=s|) 

2222  22  22  22 
Pr  ( f  j=sf/ f  j_1=Sj_1,  . . .  ,  q=s^)  Pr  (f  I^-s? r  .  .  . , 

1  2  .  * 

{ f j ; 3 >1 }  and  { f ^ ^ 1 >  Markov  chains 

11  1  1  222  2 

11  11  22  29 

pr(f  j_l=Sj_l . Pr(ff_1»sf_1, . . .  ,q=sj)  = 

1/2 

independence  of  {fjH>l}  and  {  f ^  ;  j ^.1 ) 


=Pr(-j-l=-j-l'  *  *  *  '-1=-1}  Pr C£j=s j  , _i=s j_i ) /pr  ( jf  j_i=s^_i)  (A.  5) 

The  series  of  steps  in  (A. 5)  prove  (A. 4),  The  stationarity  of 


{fj,*j>l}  follows  from  the  fact  that 


-2 


Pr(f.=sj/fj_1=sj_1)  =  (q-l)  ; j  >2  (A. 6) 
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(A. 5)  and  (A. 6)  prove  the  Lemma.  We  will  now  prove 


proposition  1. 

Let  us  first  denote  by,  T1(M),  the  following  statement. 

"The  probability,  Pr (S (1,M) =i/f|=s| , f ^=s^) ,  with  s^=s^,  is 

1  2 

independent  of  the  actual  values  of  s1  and  s^,  for  every  l  such 
that  Osi^M." 

and  by,  T2 (M) ,  the  following  statement. 

"The  probability,  Pr (S (1,M) =i/f J=s|, f i=s^) ,  with  s^#s*,  is 

1  2 

independent  of  the  actual  values  of  s1  and  s^  for  for  every  l 
such  that  O^isM." 

We  will  prove,  by  induction,  that  T1(M)  and  T2 (M)  are  true 
For  M=1  we  have: 

(A. 7) 

(A. 8) 


see  figure  1 


Pr(S(l,l)=0/fJ=sJ,f*=s*)=0 

12 

?S1  S1 

Pr(S(l,l)=l/f|=sJ,f*=s*)=l 


;s1=s1 

Pr  (S  ( 1 , 1)  =0/^=3 1 ,  f *=s*) = 

;sl *S1 

=Pr ( f {*£ l , f \*£ \/t i-. I , f i=s l ) 


'i 

E  1  21/21  1  2222  221  122 

Pr(f1#f1,f1^f2/f1=s1,f1=s1,f2-s2)  Pr(f2=s2/f1=s1,f1=s1)  = 


;s2ysl'  s2/s 
=(q-2)/(q-l) 


1 

1 


(A. 9) 
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i 


(A. 10) 


Similarly,  we  can  show  that 
Pr(S(l,l)=l/fJ=sJ,f2=s2)=l/(q-i)  ;s|fs2 

(A. 7) , (A. 8) , (A. 9)  and  (A. 10)  indicate  that  Tl(l)  and  T2(l) 
are  true.  Suppose  that  T1(M)  and  T2 (M)  are  also  true  for  M=n-l.Let 
us  now  denote 

b1(i)=Pr(S(l,n-l)=i/f|=s|/f^=s^)  ;0<i<n-l,  sj=s2  (A. 11)  and 
b2 (i)=Pr(S(l,n-l)=i/f|=s|, f2=s2)  ;0<i<n-l,  s^s2  (A. 12) 

We  will  now  show  that  T1(M)  and  T2 (M)  are  true  for  M=n.We 
first  note  that 

Pr  (S  ( 1 ,  n)  =0/f^=s11 ,  f2=s2)  =0  (A.  13  ) 

12 

;S1=S1 

Furthermore, 

Pr(S(l,n)=i/f|=s11,f^=s^)  = 

1  2  „  . 

;s1=s1,  0<i<n 

q  q 


,„1  1  „2  2,1  1  2  2 , 

Pr (f2=s2, f2=s2/f1=s1, f1=s1)= 


'l  'i 

E  E 


S2=1 


Pr(S(l,n)=i/fJ=s11,fJ=s^,f^=s2,f^=s^)  (q-l)"2  + 


21  .  2a  1 

s2=s2  and  s2#Sl 


Pr(S(l,n)=i/fJ=s* 


.  2  1 
and  s  =s, 


(q-l)'2+ 


+  ^  2_,  Pr  (S  (1 ,  n)  =i/f  1=s1 ,  f*=Sl,  ^2=s2  ’  f  2  S2)  ^-1)  + 

■4  M 


S2=1  S2=1 

2  1  2  1 

s2?<s2  and  s2^si 


4  q 

y:  y  Pr(S(l,n)=i/fJ=sJ,f*=s*,f2=s2,f*=s*)  (q-l)~2= 


Sl=1  S2=1 

s^s2  and  S2=S1 


12  2  1 

(we  cannot  have  s1=s1  and  s2=s1) 


y  y  Pr(S(2,n)=i-l/f11=s11,f^=s^,f2=s2,f^=s^)  (q-l)"2  + 


a1-!  s2=1 

s2  J-  s2  l 

2  1  ,  2/  1 

S2*S2  and  S2^S1 


q  q 

y  y  Pr (S (2 , n) =i-l/f^=s| , f^=s^, f 2=s2 , f 2=s2 ,  (q-l)"2= 

1  2 

c  =1  c?  —  1  /  ^  •  -i  /  t  ~  — _ 


s2=1  s2-l 

2  ,  1  .  2  ,  1 

s2*s2  and  s2^S;L 


is  a  Markov  chain  (Lemma  1) 


■£  £ 

a1-!  s2-! 

S2  1  S2_1 


Pr (S (2 , n) =i-l/f 2=s2 , f \=s\)  (q-1) '2+ 


_2  1 


2  ,  1 


=s2  and  s2/Sl 


'I  M 

£  £ 

■i-i  s?=i 


Pr(S(2,n)-i-l/f’=s’,f|-s*>  (q-1) 


is  stationary  (Lemma  1) 


LV.V.W^  V  .VJVW.V.  7VV  "V'-’V"  %ivvi^n7v 


^  ^  PrfSd.n-D^-l/f^s’,^)  (q-l)"2+ 


S2=1 


-;-1 


S2~S2  and  S2^S1 


2  /  I 


^  ^  Pr(S(l,n-l)=i-l/fJ=sJ,f2=s*)  (q-l)“2 


(A. 14) 


S2=1 


s^l 


and  s^sj 


We  now  make  the  following  observations. 


12  12 

0.1)  If  f1  is  in  state  (s-l/S-l)  with  s1=s1,  then  f2  can  be  in  q-1 

12  2  12  1 

different  states  (s2,s2)  with  the  properties  s2=s2  and  s^s.^, 


1  2 

independently  of  the  actual  values  of  s^  and  s^. 


12  12 

0.2)  If  f1  is  in  state  (s^s^  with  s1=s1,  then  12  can  be  in 


2  1 

(q-1) (q-2)  different  states  with  the  properties  s2^s2  and 


2  1.  12 

®2^s1 ,  independently  of  the  actual  values  of  s^^  and  s^ 


From  (A. 14),  0.1  and  0.2  we  conclude  that 


Pr(S(l,n)=i/fJ=s^,f^=s^)=b1(i-l)  (q-1) /(q-1) 2+ 

1  2  . 

;s1=s1,  ocisn 

+b2(i-l)  (q-1) (q-2)/(q-l)2 


(A. 15) 


Similarly,  we  can  show  that 


Pr(S(l,n)=0/fJ=s^,f2=s*)=(q-2)2/(q-l)2  b2(0) 


(A. 16) 


Pr(S(l,n)=i/f|=s^, f 2=s^)=  (q-2)/(q-l)2  b^i)* 

1  ,  2  . 

;s1fs1,  0<i<n 


+(q-l)/(q-l)2  b2(i-l)+  (q-2 ) 2/ (q-1) 2  b2(i) 


(A. 17) 


m 

v  *  w  * 


and, 

Pr(S(l,n)=n/fJ=sJ,f*=s*)=(q-l)/(q-l)2  b2(n-l) 

1  ,  2 

;S1*S1 


(A. 18) 


Expressions  (A. 15)  through  (A. 18)  prove  that  T1(M)  and  T2  (M) 
are  true  for  M=n  as  well.  Hence,  by  induction,  T1(M)  and  T2 (M)  are 
true  for  every  M£l. 

Furthermore,  the  induction  procedure  above  and  in  particular 
expressions  (A. 15)  trough  (A. 18)  verified  the  validity  of 
equations  (5)  through  (7)  of  section  3. 
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Table  1 


(31,7)  RS  code 


q 

Pe<2> 

Pe(2) 

IPe(2)-pe(2)l  /  P  e  (  2  ) 

10 

0.78507246 

0.77127123 

0 . 17  894 13  3D-0 1 

25 

0.93451075D-01 

0. 97 636106D-01 

0.42863559D-01 

50 

0. 66080673D-02 

0 . 72908919D-02 

0. 93654467D-01 

(31,15)  RS  code 


q 

P«<2> 

Pe(2) 

IPe  (2)  -pre  (2)1  /pe(2) 

10 

0.13915048 

0. 15076248 

0.77021816D-01 

25 

0. 36450112D-03 

0. 53081917D-03 

0.31332336 

50 

0. 14 188960D-05 

0. 23669732D-05 

0.40054412 

(31,23)  RS  code 


q 

Pe(2) 

Pe(2) 

|Pe(2)-pe(2)|  /Pe(2) 

10 

0 . 25728328D-02 

0.43874633D-02 

0.41359445 

25 

0. 9370620D-07 

0 . 28409390D-06 

0.67015764 

Table  1  (continued 


(63,15)  RS  code 


q 

Pe(2) 

Pe(2> 

|Pe(2)-frG(2)|  /pe(2) 

10 

0.91745546 

0.90631180 

0. 12295613D-01 

25 

0. 58637826D-01 

0. 62530221D-01 

0. 62248220D-01 

50 

0 . 74  3  9  5697D-03 

0. 87602105D-03 

0. 15075445 

(63,31)  RS  code 


q 

Pe(2) 

Pp(2) 

lPe(2)-?e(2)|  /pe(2) 

10 

0.10016353 

0.11197709 

0. 10549979 

25 

0. 32382994D-05 

0. 62966015D-05 

0.48570678 

50 

0. 11704D-09 

0. 29755D-09 

0.60665434 

(63,47)  RS  code 


q 

Pe(2) 

|pe(2)-pe(2)|  /pe(2) 

10 

0. 1054 14  8  5D-0  3 

0.26357027D-03 

0.6000503 

(31,15)  RS  code 


q 

Pe(3) 

Fe(3) 

ipe(3)-gre(3)|  /Pe<3) 

10 

0.85167393 

0.84020289 

0. 13652702D-01 

25 

0. 35959153D-01 

0. 39002482D-01 

0. 78029111D-01 

50 

0. 38309009D-03 

0.45761967D-03 

0.16286358 

(31,23)  RS  code 


Pe(3) 

^e(3) 

IPe(3)“Pe(3)l  /Pe(3) 

0 . 29646470 

0. 30376087 

0.24019453D-01 

0 . 2  5272  7  62  D-03 

0. 35074331D-03 

0.27945134 

0 . 13704329D-06 

0. 22482071D-06 

0.39043298 

(63,31)  RS  code 


q 

Pe(3) 

Pe(3) 

10 

0.95754921 

0.95027964 

0. 76499271D-01 

25 

0 . 10935348  D-0 1 

0. 12442345D-01 

0.12111840 

50 

0. 35134889D-05 

0 . 4809554  8D-05 

0.26947731 

q 

Pe(3) 

Pe(3) 

|pe(3)-f>e(3)|  /pe(3) 

10 

0.30594667 

0.31322755 

0. 23244998D-01 

25 

0. 149 1614  8D-05 

0. 26688242D-05 

0.44109664 

50 

0. 115D-11 

0. 286D-11 

0.59790209 
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